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ABSTRACT 

Recently, a number of ultracool dwarfs have been found to produce periodic radio bursts with high 
brightness temperature and polarization degree; the emission properties are similar to the auroral 
radio emissions of the magnetized planets of the Solar System. We simulate the dynamic spectra of 
radio emission from ultracool dwarfs. The emission is assumed to be generated due to the electron- 
cyclotron maser instability. We consider two source models: the emission caused by interaction with 
a satellite and the emission from a narrow sector of active longitudes; the stellar magnetic held is 
modeled by a tilted dipole. We have found that for the dwarf TVLM 513-46546, the model of the 
satellite-induced emission is inconsistent with the observations. On the other hand, the model of 
emission from an active sector is able to reproduce qualitatively the main features of the radio light 
curves of this dwarf; the magnetic dipole seems to be highly tilted (by about 60°) with respect to the 
rotation axis. 

Subject headings: radio continuum: stars - stars: low-mass - stars: brown dwarfs - stars: magnetic 
held - planet-star interactions 



1. INTRODUCTION 

Very low-mass stars and brown dwarfs (collectively 
known as ultracool dwarfs) represent a distinctive class 
of radio sources. Firstly, they violate the empirical 
Gude l-Benz relation between the rad io and X-ray emis- 
sion ([Berger et al.| [2005: Bcrgcr 2006). For a wide range 
of objects (including stars of various spectral classes 
and solar and stellar flares), the spectral radio lumi- 
nosity Lji and X-ray lumi nosity Lx satisfy th e rela- 
tion L X /L R ~ 10 15 - 5 ± 5 Hz (jGiidel fc Benz|[l99l which 
indicates the common origin of these emissions. For 
the stars with spectral class M7 and cooler, the X- 
ray luminosity abruptly decreases dBerger et al.l I2010T) 
as a result of transition fr om hot ionized coronae to 
cooler neutral atmospheres (Schriiverl l2009D . Neverthe- 
less, it was recently discovered that a number of very 
low-mass stars and brown dwarfs in the spectral range 
M7-L3 are the sources of unexpect edly intense radio 
emissi o n in t h e GHz freque n cy range dBerger et aLll200lL 
2008al [20101: iBergerl 120021 l2006t iBurgasser fc Putmanl 



2005UC)sten et all 120061 720091: jPhan-Bao et al.l 120071; 
Antonova et al.ll2008HRavi et al.ll201lD ; the emission in- 



tensity is comparable to or even exceeds that of the early 
M dwarfs. Secondly, the radio emission from ultracool 
dwarfs can include (in addition to a quiescent compo- 
nent) short periodic pulses with high brightness temper- 
ature and almost 100% circular polarization, whose pe- 
riod seems to coinc i de wit h the s tellar rotation period 
(IBerger et all 120051. l27)08bl [20091 : IHallinan et all l200l 
2007L 120081 iDovle et al.1 120101 ) . The total pulse dura- 
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tion can be as short as At < 10 2 T ro t (Hallina n et al.1 
l2007tlBerger et al.ll2008bll2009f ). where T rot is the stellar 
rotation period, which suggests high emission directivity 
and small (relative to the stellar radius) source size. The 
emission pulses can exhibit fine temporal structure that 
is repeated each rotation period, with the duration of in- 
divid ual subpulses down to a few seconds (Hallinan et al. 
[2001 . 

The above features are untypical for the stellar radio 
emission and more similar to the auroral radio e mission 
of th e magnetized planets of the Solar System (jZarkal 
1998), but with much higher frequency and intensity. 
In particular, the most likely emissio n mechanism is 
the electron-cyclotron maser instability ()Wu fc Ledll979l : 
iMelrose fc Dulklll982t rTreumannl 120061) which produces 
emission near the local electron cyclotron frequency. In 
order to provide the observed emission frequencies, ul- 
tracool dwarfs must possess magnetic fields of several 
thousand gauss. Observations of Zeeman broadening of 
molecular lines in cool stars (with the spectral class up 
to M9) have confirmed the existen ce of surface magnetic 
fields with strengths up to 3900 G (jReiners fc Basriil2007l 
2010; IReiners et all [20091) . 

The interiors o f ultracool dwarfs are ex pected to be 
fully convective (jChabrier fc Baraffj 0,9971 ) . and there- 
fore the magnetic dynamo processes in these objects 
should differ significantly from those in the Sun and other 
solar- type stars; seyeral models have be e n pro p osed (e.g., 
Dobler et all 120061 : iChabrier fc Kiikerl 120061: iBrowningl 



2008). The observations have shown that the morphol- 



ogy of the periodic light curves and hence the magnetic 
field stru cture remain stabl e at timescales of up to several 
months ([Doyle et al.ll2010l ). although observations sepa- 
rated by longer times (years) can dem onstrate consider- 
able variations ([Antonova et al.l 120071 ). Such behaviour 
can be reconciled with a turbulent dynamo model. 

The magnetic field topologies of the radio-emitting ul- 
tracool dwarfs are unknown. The planets of the So- 
lar System have, in general, dipole-like magnetic fields 
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Figure 1. Sketch of the emission source models: a) the model 
with the star-satellite interaction; b) the model with emission from 
a set of magnetic field lines (the rotational and magnetic equatorial 
planes are shown by dotted lines). The magnetic dipole axis is 
shown by bold arrows. 

with a dipole either nearly parallel to the rotation 
axis (as for the Earth, Jupiter, and Saturn) or highly 
tilted with res pect to the rotation axis (a s for Uranus 
and Neptune) (jRussell fc Dougherty! 1201(3) . Spectropo- 
larimetric observa tions of a sample of M5-M8 dwarfs 
(M orin et aT]|2010D have revealed that these stars tend to 
have dipolc-like axisymmetric magnetic fields (although 
other topologies have been observed as well). On the 
other hand, s imulations of the a 2 dynamo in fully con- 
nective stars (IKiiker fc Rudiger| [l99"a IDobler et~aUl2006t 
IChabrier fc Kiikerl 120061) predict approximately dipole- 
like magnetic field with the dipole perpendicular to the 
rotation axis. Observations of the cyclotron-maser ra- 
dio emission offer a unique (and often the only) way to 
study the magnetic field topology and other parameters 
of the magnetospheres of ultracool dwarfs; the obtained 
results, in turn, will be very important for the further 
development of the stellar dynamo theory. 

The aim of this work is to analyze and compare dif- 
ferent source models that can result in formation of the 
short periodic (or quasi-periodic) emission pulses from 
ultracool dwarfs. Namely, we consider modulation of the 
emission by the stellar rotation (in a non-axisymmetric 
magnetosphere model) and the emission induced by a 
satellite orbiting the ultracool dwarf; both models are 
similar to the cases observed in the Solar System. The 
used models are described in Section [2j The simulation 
results are presented in Section [3] and compared with the 
observations in Section 01 The conclusions are drawn in 
Section [5] 

2. MODEL 

We have developed a numerical code to simulate the 
dynamic spectra of the radio emission. It is assumed that 
the magnetic field of the ultracool dwarf is dipole-like 
with the dipole located at the star center (no offset) and 
tilted by the angle S relative to the rotation axis. The 
rotation axis is inclined relative to the line of sight by 
the angle i. Since the periodic pulses can be very short, 
their source should occupy a narrow range of longitudes. 
We therefore consider two models that can provide the 
required source geometry: 

• Satellite-induced radio emission (analogous to the 
Io- Jupiter system). In this model, a satellite (ex- 
oplanet) orbits the ultracool dwarf and its inter- 
action with the magnetospheric plasma results in 
electron acceleration (possibly, due to the excita- 



tion of Alfven waves). Thus we assume that the ra- 
dio emission is generated at the satellite flux tube, 
i.e., at the magnetic field line passing through the 
satellite (see Fig. QJi). Since the Alfven waves 
propagate at a finite speed, rapid rotation of the 
dwarf can shift the radio-emitting region ahead of 
the satellite; however, this effect is neglected in this 
study. We also assume that the satellite orbit is 
circular and the orbital plane coincides with the 
stellar equatorial plane. 

• Radio emission from a sector of active longitudes 
(similar to the case of Jovian hectometric radia- 
tion). In this model, the emission is generated 
at fixed (in the rotating frame of the star) mag- 
netic field lines. Namely, we assume that the 
"active" field lines cross the magnetic equatorial 
plane at a given distance from the star center (i.e., 
have the same L-shell numbers) and their magnetic 
longitudes A are confined within a narrow range 
A A <C 1. Figure [Tp demonstrates two sets of "ac- 
tive" lines: with A ~ 0° and A ~ 90° (the mag- 
netic longitude is calculated relative to the plane 
containing the rotation axis and the dipole axis). 
The electrons are accelerated, most likely, due to 
magnetic reconnection at the boundary of the mag- 
netospheric corotation region. 

Similar models (but with significantly diff erent parame- 
ters) were studied by FIcss & Zarkal (|2011[ ) in application 
to the exoplanetary or exoplanet-induced stellar radio 
emission. 

The emission is produced due to the electron-cyclotron 
maser instability. We consider two types of unsta- 
ble electron distributions: the loss c one and the she ll 
(or horseshoe) distribution (see, e.g.. iTreumannl [20061. 
The shell distributions have been observed in situ 
within the sources of terrestrial and Saturnian kilomet- 
ric radiation (e.g ., iDelorv et aT]|1998t lErgun et all 120001: 
lLamv et al.ll2010l ) , while the characteristics of Jovian de- 
cametric radiation agree better with the predictions for 
the loss-cone-drive n instability (jHess et al.l 120071 120081: 
IMottez eT al. 2010). Thus both types of distributions can 
occur (possibly, simultaneously) in the magnetospheres 
of ultracool dwarfs. For the loss-cone-driven instabil- 
ity, the emission frequency / and propagation direction 
9q (relative to the local magnetic field) are given by 
(jKuznetsov fc Tsadl2007l: IHess fc Zar ka 2011) 



cos 0o — i v / c )/ cosa c , 



(1) 

(2) 

where /b is the local electron cyclotron frequency, v is 
the typical speed of the energetic electrons, and a c is 
the loss-cone boundary. The emission is assumed to be 
azimuthally symmetric with respect to the local magnetic 
field, thus it is beamed along a thin conical sheet with 
the half- apex angle of 6q. The loss-cone boundary follows 
the transverse adiabatic invariant 



sin 2 a c = /b//bd 



(3) 



where /Bmax is the electron cyclotron frequency at the 
footprint of the radio-emitting magnetic field line (where 
it enters dense atmospheric layers). 
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For the shell-dri ven instability, the emission fre- 
quency is giv en _by (|Hess et al.l 120081 iKuznetsovl l20lH: 
IHess fc ZarlfflllMll) 

/~/bV1-« 2 /c 2 </b, (4) 

and the emission direction 9q ~ 90° irrespective of the 
electron energy or the source location. 

We compute the dynamic spectra in the following way: 
for a given time, frequency, and radio-emitting magnetic 
field line (the line is either fixed or determined by the cur- 
rent satellite position relative to the magnetic dipole) , we 
find the emission source coordinates. Then we calculate 
the angle 6 between the local magnetic field and the line 
of sight. Contribution of the considered field line to the 
emission intensity is calculated as 

I(tJ)~exp[-(e-6 ) 2 /A6 2 ], (5) 

where AO is the beam half-width. The radio emission 
is assumed to be 100% circularly polarized; the emis- 
sions originating from the northern and southern mag- 
netic hemispheres have opposite polarization signs. As 
Eq. ([5]) implies, the considered model does not take into 
account variations of the emission intensity with height; 
therefore it provides only some general emission patterns 
in the time-frequency domain, i.e., where and when the 
pulses can occur. 

Radio emission can propagate only if its frequency 
exceeds the local plasma frequency, which can result 
in an additional tim e-dependent low-frequency cutoff 
(|Hess fc Zarkal 1201 ID . Due to high gravity and low 
plasma temperature (as indicated by very weak X-ray 
emission), the plasma density scale height in the mag- 
netospheres of ultracool dwarfs is expected to be much 
less than the stellar radius. Therefore we consider a 
simplified occultation effect: the emission (regardless of 
frequency) is completely absorbed below some critical 
distance from the star center R cr and propagates freely 
above that distance; we assume that R cr approximately 
equals the stellar radius i?*. The other propagation ef- 
fects (such as refraction, scattering, and depolarization) 
are assumed to be negligible. 

3. SIMULATION RESULTS 

Even within the above mentioned simplified models, 
it is impossible to consider here all the possible combi- 
nations of the source parameters. Therefore we restrict 
ourselves to searching for the cases when the simulated 
radio emission would be similar to the observed one. As 
an example, we use the most studied radio-emitting ul- 
tracool dwarf TVLM 513-46546 (hereafter TVLM 513). 
T his M9 dwarf has th e rotation period of T rot = 1.96 
h (|Hallinan et alj|2007t iDovle et all 120101 ): its mass and 
radius can be estimated as M* = (0. 06 — 0.08)M(yi and 
R* = (0.097- 0.109) ^, respectively (|Dahn et alJl2002b 
Mailman et al.l I20061 ). The obs erved projected rotatio n 
velocity of v sin i ~ 60 km s" 1 (jMohantv fc Basrill2003] ). 
together with the mentioned radius and rotation period, 
imply that the rotation axis should be hig hly inclined rel- 
ative to the line of sight ( 65° < i < 90°) (jHallinan et al.l 
120061 ; iBerger et ai]|2008b[ ). In the below simulations, we 
use the following average values of the dwarf mass and 
radius: M* = O.O7M , R* ~ O.LR = 70 000 km. 



In the simulations, we assume that the radio beam half- 
width is A9 = 2° and the speed of the energetic electrons 
is v — 0.3c. The electrons responsible for the planetary 
radio emissions in the Solar System typically have lower 
energies (v ~ 0.1c), but we chose a higher value in order 
to make the differences between the loss-cone-driven and 
shell-driven instabilities more apparent. 

3.1. Satellite-induced emission 

Since we are interested in a periodic radio emission 
with the period equal or close to the stellar rotation pe- 
riod T ro t, we consider the case when the satellite moves 
relatively slowly, so that its orbital period T or b 3> T rot . 
In the below simulations, we assume that T or b = 20.1T rot ; 
for the TVLM 513 system, this corresponds to the or- 
bit radius of R OI b = 24. 0i?». If the magnetic dipole is 
axisymmetric, the emission pulse period will equal the 
satellite orbital period; therefore we consider only the 
nonzero dipole tilts. 

Figure [2] shows the simulated dynamic spectra (emis- 
sion intensity) for low (5 = 15°) and high (6 = 60°) 
dipole tilts and for two inclinations of the rotation axis 
(i = 90° and 60°). The emission frequency in the spectra 
is expressed in units of the electron cyclotron frequency 
at the magnetic pole /bo- We have found that if the mag- 
netic dipole tilt is relatively small, then the emission can 
be detected (i.e., it is beamed towards an observer) only 
during some favourable orbital phases of the satellite. As 
a result, the emission pulses at a single frequency occur in 
series separated by periods of inactivity; the favourable 
orbital phases are frequency-dependent and depend also 
on the rotation axis inclination and the characteristics of 
the electron distribution function. 

For the larger dipole tilts, wide-band emission pulses 
can be detected at any orbital phase of the satellite. The 
required tilt increases with a decrease of the rotation 
axis inclination and can be estimated as S > 120° — i 
(we consider here only the inclinations close to 90°). 
The dynamic spectra of the shell-driven emission look 
like series of nearly vertical stripes (with high frequency 
drift rates); the high-frequency cutoff is determined by 
the electron cyclotron frequency at the footprint of the 
satellite flux tube /emax with the relativistic correction 
^ . For the loss-cone-driven emission, the dynamic spec- 
tra of separate pulses have an "inverted U" shape; the 
high-frequency cutoff is determined mainly by the visi- 
bility conditions since the emission frequency and prop- 
agation direction are interrelated. At a single frequency 
we can see, in general, four emission pulses (two with 
the right circular polarization and two with the left cir- 
cular polarization) per each rotation period; the pulses 
with opposite polarizations can overlap thus forming a 
pulse with a higher intensity and lower polarization de- 
gree. If the rotation axis inclination differs from 90° then 
the pulses with opposite polarizations become asymmet- 
ric (e.g., they have different high-frequency cutoffs); this 
effect is negligible for the shell-driven emission but sig- 
nificant for the loss-cone-driven emission. 

The emission pulses are not strictly periodic; however, 
each pulse is repeated approximately one rotation period 
later. Time intervals between the corresponding pulses in 
the adjacent rotation periods are shown in Fig. [3] (these 
intervals are nearly independent on the source param- 
eters provided that the dipole tilt is sufficiently high). 
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Figure 2. Model dynamic spectra (Stokes I) of a satellite-induced radio emission for the different inclinations of the rotation/orbital axis 
i and magnetic dipole tilts 8. a) Shell-driven emission, b) Loss- cone- driven emission. 
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Figure 3. Time intervals between the emission pulses with the 
same polarization at the fixed frequency of / = 0.6/bo (for the 
model with a satellite, i = 90°, 8 = 60°, and shell-driven insta- 
bility). Only the pulses with odd numbers are considered (i.e., 
approximately one pulse per rotation period). Diamonds and tri- 
angles correspond to the pulses with positive (V > 0) and negative 
(V < 0) circular polarizations, respectively; the abscissa represents 
the centers of the time intervals. 



One can see that the time intervals are very close to the 
stellar rotation period. Variations of the dynamic spec- 
tra and light curves from period to period can fall below 
the time resolution of the instrument if the radius of the 
satellite orbit is sufficiently large. Therefore, if observa- 
tions cover only a few consecutive rotation periods, one 
may conclude that the emission is periodic. However, 
time intervals between the emission pulses are not con- 
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Figure 4. Model dynamic spectra (Stokes J and V) of a satellite- 
induced radio emission for i = 90° and 8 = 60° ; two time intervals 
corresponding to the different orbital phases of the satellite arc 
shown. The areas which are darker or brighter than the background 
correspond to the positive or negative values, respectively, a) Shell- 
driven emission, b) Loss-cone-driven emission. 



stant; moreover, the intervals between the consecutive 
left-polarized pulses vary in antiphase with the intervals 
between the right-polarized pulses. Therefore, the dy- 
namic spectra and light curves for a chosen rotation pe- 
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riod vary significantly on time scales comparable to the 
orbital period of the satellite. This can be seen in Fig. [4] 
which shows enlarged fragments of the dynamic spectra 
in intensity and polarization for the orbital phases of the 
satellite near 0° and 90° (relative to the line of sight). 

3.2. Emission from an active sector 

In this model, the emission is always periodic with the 
period equal to the stellar rotation period. The emis- 
sion is generated at the magnetic field lines connected to 
the particle acceleration region that, in turn, is located 
in the equatorial plane at (or near) the magnetospheri c 
corotation radius (|Andre et al.lll988l ; lLinskv et al.l 1992). 
Due to fast rotation, the corotation radii and hence the 
L-shell numbers of the radio-emitting field lin es for ultra- 
cool d warfs are expected to be rather small; iRavi et ahl 
(|2011| ) estimated the corotation radius for the radio- 
emitting M8 dwarf DENIS 1048-3956 as R c ~ (1.2 - 
2.9)i?*. In the below simulations, we assume that the 
"active" magnetic field lines have L = 2. We also assume 
that the angular width of the active sector is AA = 10°; 
the total emission is calculated as the sum of emissions 
from ten field lines evenly distributed within the range 
AA. 

We should note that the existence of a narrow sec- 
tor of active longitudes implies that the magnetospheric 
topology is more complicated than the dipolar model de- 
scribed in Section [2] e.g., the dipole can be offset from 
the rotation axis (like at Neptune) and/or the higher- 
order components of the magnetic field can be important. 
Nevertheless, we assume that the magnetic field is still 
approximately dipole-like, but the conditions favourable 
for the emission generation are met only at certain field 
lines. 

Figure [5] shows the simulated dynamic spectra (emis- 
sion intensity and polarization) for low (6 — 15°) and 
high (5 = 60°) dipole tilts and for two inclinations of 
the rotation axis (i = 90° and 60°). The active sector is 
located around the magnetic longitude Ao = 0°, i.e., the 
"active" field lines are parallel to the local meridian (see 
Fig. [lb)- Since the emission source is symmetric with 
respect to the plane containing the rotation axis and the 
dipole axis, the dynamic spectra are also symmetric and 
the number of emission pulses (at a given frequency) per 
rotation period is always even. For large rotation axis 
inclinations (i ~ 90°) and low dipole tilts (S < 15°), 
the emission sources located in opposite magnetic hemi- 
spheres are observed nearly simultaneously which results 
in weakly-polarized emission. With an increase in the 
dipole tilt, the emission pulses produced in opposite 
hemispheres become asymmetric and separated in fre- 
quency (so that only the left- or right-polarized emission 
is observed in certain frequency ranges). A decrease in 
the rotation axis inclination has a similar effect, but also 
introduces a separation in time. The loss cone produces 
the "inverted-U" -shaped and ring-shaped structures in 
the dynamic spectrum, while the shell distribution pro- 
duces either vertical stripes (for i = 90°) or the pulses 
with a frequency drift (for other inclinations). For high 
dipole tilts and small rotation axis inclinations, only the 
emission from one of the magnetic hemispheres can be 
observed (this effect is more pronounced for the loss- 
cone-driven emission). 

Figure [6] shows the simulated dynamic spectra for the 



case when the active sector is located around the mag- 
netic longitude A = 90° (see Fig. [Tp). The magnetic 
field lines with A = 90° (or 270°) cross the stellar equa- 
torial plane at the maximal distance from the rotation 
axis (among the lines with a given L), which is expected 
to result in the most efficient particle acceleration. Like 
in the previous case, the loss cone produces the "inverted- 
U" -shaped and ring-shaped pulses that exhibit a reversal 
of the frequency drift at the maximal frequency, while the 
shell distribution produces pulses with a simpler struc- 
ture, whose high-frequency cutoff is determined by the 
maximal possible emission frequency at the "active" field 
line and by the occultation effects. At a given frequency, 
we can observe up to four pulses (two with the right cir- 
cular polarization and two with the left circular polariza- 
tion) per rotation period. For the loss-cone-driven emis- 
sion, the number of emission pulses at a given frequency 
is always even; for the shell-driven emission, the occulta- 
tion effects can remove some pulses thus resulting in one 
or three pulses per rotation period. Such asymmetric 
light curves can occur if the rotation axis is not perpen- 
dicular to the line of sight (i < 90°). Both an increase in 
the dipole tilt and a decrease in the rotation axis inclina- 
tion result in an increasing time separation between the 
left- and right-polarized pulses, but these factors have 
different effects on the light curves; e.g., if 6 > (90° — i), 
the emission pulses at a single frequency have a "left- 
right-left-right" pattern, while for 8 < (90° — i) we will 
obtain a "left -right -right-left" pattern. 

4. COMPARISON WITH OBSERVATIONS 

Most observations of the radio emission from ultracool 
dwarfs have been performed at one or two frequencies. 
Only a few dynamic spectra (with a rather narrow fre- 
quency range) have been recorded to date. Therefore 
we focus here on the characteristics of the emission light 
curves at a single frequency. 

TVLM 513 has the longest total duration of radio ob- 
servations among all ultracool dwarfs, including continu- 
ous observation sess ions covering several consequent ro- 
tation periods. E.g.. lHallinan etaD |2007) presented the 
observations that lasted for ~ 10 h or about five rota- 
tion periods (see Fig. 1 in the cited work for the light 
curves); the conclusion about a stable period was made. 
We have re-examined those observations and found that 
the relative deviation of a pulse period from the average 
value does not exceed 5 x 10~ 3 (for both polarizations). 
Even more strict limit ation on the possib le period varia- 
tions was obtained by iDovle et al.l (|2010f) , who analyzed 
two data sets (each lasting for ~ 8 h) separated by ~ 40 
days; the period was found to be stable with an accuracy 
of up to 10~ 5 . These results cannot be reconciled with 
the model of satellite-induced emission unless the satel- 
lite orbital motion is very slow. However, according to 
our simulations, to obtain the relative variations of the 
pulse periods of < 5 x 10~ 3 one needs the satellite orbital 
period of T or b/T rot > 100 and thus the orbit radius of 
Roib/R* ^ 70. It is highly unlikely that so distant satel- 
lite would be able to modulate the stellar radio emission, 
and therefore the model of satellite-induced emission can 
be ruled out for TVLM 513. Nevertheless, this model can 
work at other ultracool dwarfs. 

It is interesting to note that iForbrich fc Bergerl {2009) 
suggested the detection of a companion of TVLM 513 at 
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Figure 5. Model dynamic spectra (Stokes / and V) of the radio emission from an active sector for the different inclinations of the rotation 
axis i and magnetic dipole tilts 5. The areas which are darker or brighter than the background correspond to the positive or negative 
values, respectively. The "active" field lines have the magnetic longitudes around A ~ 0°. a) Shell-driven emission, b) Loss-cone-driven 
emission. 



a distance of about 20i?* from the star (i.e., with the pa- 
rameters close to those considered in Section l3Tj) . How- 
ever, it seems that this companion does not affect the 
periodic cyclotron-maser radio emission from the dwarf 
(at least, it does not affect the pulse occurrence times). 
Most likely, this is because the satellite, if present, is 
located beyond the corotation radius of the stellar mag- 
netic field. 

Figure [7] presents an example of time profiles of the 
radio emission from TVLM 513 at different frequencies; 
due to techni cal reasons, onl y the Stokes V radio flux 
is shown (see Hallinan 120091 for a more complete dis- 
cussion of the observational data). The emission was 
observed with the Wide-band Arecibo Pulsar Processor 
at the Arecibo Observ atorJI] on 2008 May 19 (see also 
lAntonova et al.l[2010al| bl). The time interval shown cor- 
responds to one stellar rotation period (1.96 h). One 
can see three emission pulses (A, B, and C), with the 
pulses A and C being positively polarized (V > 0) and 
the pulse B being negatively polarized (V < 0). Time 
intervals between the pulses are — *a — 0.075T rot and 
tc — t& — 0.48T rot , i.e., it is likely that the pulses A and 
C are produced in the same source due to the shell-driven 
instability. The same morphology of the light curves was 

7 The Arecibo Observatory is part of the National Astronomy 
and Ionosphere Center, which is operated by Cornell University 
under a cooperative agreement with the National Science Founda- 
tion. 



observed in each rotation period over three consecutive 
nights (2008 May 18, 19, and 20), although the pulse 
amplitudes and durations could vary significantly (e.g., 
the pulse C in Fig. [7] is much weaker than the pulses A 
and B, but an opposite relation was observed in some pe- 
riods). A stable temporal structure with a pronounced 
periodicity implies that the emission source is fixed in 
the rotating frame of the star. 

The facts that there is an odd number of emission 
pulses per rotation period and the time interval between 
the pulses with the same polarization sign is close to 
T rot /2 suggest that (i) the emission comes from an ac- 
tive sector rotated with respect to the plane containing 
the rotation axis and the dipole axis (e.g., with A ~ 90° 
or 270°), (ii) the emission is produced due to the shell- 
driven instability, and (iii) the rotation axis inclination i 
differs from 90°. As a result, the counterpart of the pulse 
B is invisible because the emission source is occulted by 
the star at that time. 

Figure [8] shows the simulated dynamic spectrum and 
light curve for the following simulation parameters: mag- 
netic field strength at the magnetic pole Bo = 2560 G 
(which corresponds to the electron cyclotron frequency of 
/bo = 7.2 GHz), rotation axis inclination i = 70°, mag- 
netic dipole tilt S = 60°, "active" magnetic line number 
L = 2.15, active sector longitude Ao = 270°, active sec- 
tor width AA = 1°, and speed of the energetic electrons 
(with the shell distribution) v = 0.1c. The parameters 
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Figure 6. Same as in Fig. \E\ for the "active" field lines with the magnetic longitudes around A ~ 90° 
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Figure 7. Time profiles of the radio emission (Stokes V) from 
TVLM 513 observed on 2008 May 19. The intensity scales in the 
first and second halves of the considered time interval are differ- 
ent in order to make the pulse C more visible. The gaps in the 
time profiles correspond to the time intervals when the instrument 
calibration was performed. 



are chosen so as to provide an agreement with the obser- 
vations. Namely, our aim was to reproduce a three-pulse 
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Figure 8. Top: model dynamic spectra (Stokes V) of the radio 
emission from an active sector. Bottom: the corresponding light 
curve at the frequency of 4.6 GHz. The sumulation parameters are 
given in the text; they were chosen so as to provide a qualitative 
agreement of the simulated light curve with the observations. 



structure in a sufficiently wide frequency range as well as 
the above mentioned time intervals between the pulses. 
Due to the large number of free parameters, the solu- 
tion presented in Fig. |8]is not unique. However, we have 
found that only relatively high dipole tilts (S ~ 50° -70°) 
can be reconciled with the observations. Also, we have 
estimated the favorable range of the "active" magnetic 
line numbers as L ~ 2.0 — 2.6, which corresponds to the 
surface magnetic field values of B ~ 2800 — 2100 G. 
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Despite a qualitative agreement, there are apparent 
differences between the simulations and the observations. 
Firstly, the observed pulses exhibit much higher fre- 
quency drift rates than predicted by the simulations. 
Possibly, this is because either the magnetic field is not 
exactly dipole-like or the emission sources do not fol- 
low the magnetic field lines (e.g., due to the peculiarities 
of the particle acceleration and precipitation processes) . 
Secondly, the pulses often exhibit fine spectral struc- 
ture that varies from period to period. Such fine struc- 
tures are typical of the planetary radio emissions and are 
caused, most likely, by small-scale irregular plasma inho- 
mogeneities (turbulence) in the radio emission sources. 
Therefore, dynamic spectra covering both wide frequency 
range and long time intervals are necessary to recover re- 
liably the magnetic field topology. 

It is interesting to note that the pulse C in Fig. [7] 
has a much lower intensity but a much longer duration 
than pulse A, so that the total radiated energy is ap- 
proximately the same. If we assume that these pulses 
are produced in the same source then the variations of 
the pulse intensity and duration can be caused by the 
radiation scattering during propagation, which is differ- 
ent for the different source orientations. An alternative 
explanation is that the active sector width AA can vary 
significantly over timescales of about one hour, while the 
total number of the energetic electrons in this sector re- 
mains nearly constant. 

5. CONCLUSION 

We have simulated the dynamic spectra of radio emis- 
sion from ultracool dwarfs. The stellar magnetic field 
was modeled by a tilted dipole. Two source models were 
considered: the model in which the emission was pro- 
duced as a result of the star-satellite interaction, and the 
model with emission from a narrow sector of magnetic 
longitudes. The emission was assumed to be generated 
due to the electron-cyclotron maser instability; two un- 
stable electron distributions (the loss cone and the shell 
distribution) were considered. We have compared the 
simulation results with the observations of the ultracool 
dwarf TVLM 513. We have found that: 

• The model of the satellite-induced emission is in- 
consistent with the observations since it predicts 
variable time intervals between the emission pulses, 
whereas the radio emission from TVLM 513 ex- 
hibits a stable period. 

• The model of emission from an active sector is able 
to reproduce the main features of the radio light 
curves of TVLM 513 (such as the pulse polariza- 
tions and occurrence times). It seems that the 
magnetic dipole is highly tilted (by about 60°) with 
respect to the rotation axis, i.e., the magnetic field 
topology of the dwarf is similar to that of Uranus. 
Also, the model with the shell-like (or horseshoe- 
like) electron distribution fits the observations bet- 
ter than the model with the loss cone. 

• On the other hand, the considered model of emis- 
sion from an active sector seems to be oversimpli- 
fied since it cannot account for, e.g., the observed 
frequency drift rates, intensity variations with fre- 



quency, and variations of the pulse intensity and 
duration from period to period. 



The experimental and theoretical studies of radio emis- 
sion from ultracool dwarfs are only at their beginning. 
Dynamic spectra of the emission (with high spectral and 
temporal resolution) covering wide frequency range and 
long time intervals are necessary to make more definite 
conclusions about the magnetic field topology, particle 
acceleration mechanism, and other parameters of the 
stellar magnetospheres. 
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